%FFT频谱分析

%离散序列分析：128点DFT
clear;
N=128;
n=0:N-1;
xn=0.9*sin(2*pi/N*n)+0.6*sin(6*pi/N*n);
XK=fft(xn,N);
%求幅度谱
magXK=abs(XK);
%求相位谱
phXK=angle(XK);
subplot(1,2,1);

%画时域信号
plot(n,xn);
xlabel('n');
ylabel('x(n)');
title('x(n) N=128');
subplot(1,2,2);
%画频谱
k=0:length(XK)-1;
stem(k,magXK,'.');
% disp(XK);
xlabel('k');
ylabel('X(K)');
title('X(K) N=128');

[H,W] = freqz(xn, 1, 512);
plot(W/pi, 20*log10(abs(H)/abs(H(1))));


% 连续信号采样分析
% clear;
% %时域采样周期
% T = 1.0/3000;
% %采样点数
% N = 120;
% n = 0:T:(N-1)*T;
% k = 0:1:N-1;
% %时域信号
% xn = (1+cos(2*pi*100*n)).*cos(2*pi*600*n);
% XK = fft(xn, N);
% %幅度谱
% magXK = abs(XK);
% subplot(2,1,1);
% %对应的数字频率
% f = 2*k/N;
% stem(f, magXK);
% xlabel('\omega/\pi');
% ylabel('|X(k)|');
% title('N=120');

